New effective treatment of the light-front nonvalence contribution in 
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We discuss a necessary nonvalence contribution in timelike exclusive processes. Fol- 



lowing a Schwinger- Dyson type of approach, we relate the nonvalence contribution to 
an ordinary light-front wave function that has been extensively tested in the spacelike 
exclusive processes. A complicate four-body energy denominator is exactly cancelled 
in summing the light-front time-ordered amplitudes. Applying our method to K13 and 
D° — > K~l + vi where a rather substantial nonvalence contribution is expected, we find not 
only an improvement in comparing with the experimental data but also a covariance(i.e. 
frame-independence) of existing light-front constituent quark model. 

With the wealth of new and upgraded B-meson factories, exclusive decay processes 
will be studied intensively. Unlike the leading twist structure functions measured in deep 
inelastic scattering, such exclusive channels are sensitive to the structure of the hadrons 
at the amplitude level and to the coherence between the contributions of the various 
quark currents and multi-parton amplitudes. The central unknown required for reliable 
calculations of weak decay amplitudes are thus the hadronic matrix elements. 

Perhaps, one of the most popular formulations for the analysis of exclusive processes 
involving hadrons may be provided in the framework of light-front (LF) quantization [|TJ . 
In particular, the Drell-Yan-West (q + = q° + q 3 = 0) frame has been extensively used 
in the calculation of various electroweak form factors and decay processes |2|-|4]]. As an 
example, only the parton-number-conserving (valence) Fock state contribution is needed 
in q + = frame when the "good" component of the current, J + or J^ = (J x , J y ), is used 
for the spacelike electromagnetic form factor calculation of pseudoscalar mesons. The LF 
approach may also provide a bridge between the two fundamentally different pictures of 
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Figure 1. Effective treatment of the LF nonvalence amplitude. 



hadronic matter, i.e. the constituent quark model (CQM) (or the quark parton model) 
closely related to the experimental observations and the quantum chromodynamics (QCD) 
based on a covariant non-abelian quantum field theory. The key of possible connection 
between the two pictures is the rational energy-momentum dispersion relation (i.e. the 
square root operator does not appear) that leads to a relatively simple vacuum structure. 
There is no spontaneous creation of massive fermions in the LF quantized vacuum. Thus, 
one can immediately obtain a constituent-type picture, in which all partons in a hadronic 
state are connected directly to the hadron instead of being simply disconnected excita- 
tions (or vacuum fluctuations) in a complicated medium. A possible realization of chiral 
symmetry breaking in the LF vacuum has also been discussed in the literature ||. 

On the other hand, the analysis of timelike exclusive processes has remained as a 
rather significant challenge in the LF approach. In principle, the q + 7^ frame can 
be used to compute the timelike processes but then it is inevitable to encounter the 
particle-number-nonconserving Fock state (or nonvalence) contribution. The main source 
of difficulty in CQM phenomenology is the lack of information on the non- wave-function 
vertex(black blob in Fig. 1(a)) in the nonvalence diagram arising from the quark- ant iquark 
pair creation/annihilation. The non-wave-function vertex(black blob) was recently also 
called the embedded state |J. This should contrast with the white blob representing 
the usual LF valence wave function. In principle, there is a systematic program laid out 
by Brodsky and Hwang |7J] to include the particle-number-nonconserving amplitude to 
take into account the nonvalence contributions. However, the program requires to find 
all the higher Fock-state wave functions while there has been relatively little progress in 



computing the basic wave functions of hadrons from first principles. Recently, a method 
of analytic continuation from the spacelike region has also been suggested to generate 
necessary informations in the timelike region without encountering a direct calculation of 
the nonvalence contribution ||. Even though some explicit example has been presented 
for manifestly covariant theoretical models, this method has not yet been implemented to 
more realistic phenomenological models. 

In this letter, we thus present an alternative way of handling the nonvalence contri- 
bution. Our aim of new treatment is to make the program more suitable for the CQM 
phenomenology specific to the low momentum transfer processes. Incidentally, the light- 
to-light (Ki 3 ) and heavy-to-light (D° — > K~t r V() decays involving rather low momentum 
transfers bear a substantial contribution from the nonvalence part and their experimental 
data are better known than other semileptonic processes with large momentum transfers. 
Including the nonvalence contribution, our results on Kg 3 and D° — > K~£ + i>g not only 
show a definite improvement in comparison with experimental data but also exhibit a 
covariance (i.e frame-independence) of our approach. 

The semileptonic decay of Qiq bound state with four-momentum Pf and mass Mi into 
another Q 2 q bound state with P 2 and M 2 is governed by the weak current, viz., 

J"(0) = (P2\Q 2 YQi\Pi) = U(q 2 )(Pi + P2Y + /-(aV, (1) 

where q^ = (Pi — P^ is the four-momentum transfer to the lepton pair (lv) and mj < 
q 2 < (Mi — M 2 ) 2 . The covariant three-point Bethe-Salpeter (BS) amplitude of the total 
current J M (0) in Eq. ([[]) may be given by 

r d 4 k H™ V H™ V S» 1 

( ) - lN cJ ( 27r )4 (p2 _ m 2 + %e) ( p 2 _ m 2 + %e) (p2 _ m 2 + ^ > ( 2 ) 

where N c is the color factor, H^?Z is the covariant intitial [final] state meson-quark vertex 
function that satisfies the BS equation, and S^ = Tr[7 5 (^ 1 +mi)7 M (^ 2 -|-m2)75(— rfq+mq)]- 
The quark momentum variables are given by pi — Pi — k, P2 = P2 — k, and Pq = k. As 
shown in the literature ||, the LF energy integration reveals an explicit correspondence 
between the sum of LF time-ordered amplitudes and the original covariant amplitude. 
For instance, performing the k~ pole integration, we obtain the LF currents, Jy and J^ v 
corresponding to the usual LF valence diagram and the nonvalence diagram shown in 
Fig. 1(a), respectively. Since H™ satisfies the BS equation, we iterate H%™ once and 
perform its LF energy integration to find the corresponding LF time-ordered diagrams 
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Figure 2. Non-wave-function vertex(black blob) linked to an ordinary LF wave func- 
tion(white blob). 



Figs. 1(b) and 1(c) after the iteration. The similar idea of iteration in a Schwinger-Dyson 
(SD) type of approach was presented in Ref. |J to pin down the LF bound-state equation 
starting from the covariant BS equation. Comparing the LF time-ordered expansions 
before and after the iteration, we realize that the following link between the non-wave- 
function vertex (black blob) and the ordinary LF wave function (white blob) as shown in 
Fig. 2 naturally arises, i.e., 

(M 2 - M^'( Xi ,k ±i ) = y[dj/][d 2 l±]/C(a; i ,k J . i ;j/ i ,l ij -)*(j/ J -,l ij -), (3) 

where M is the mass of outgoing meson and Mq = (ml + k^J/xi — (m 2 , + k 2 L2 )/(— x 2 ) 
with X\ = 1 — X2 > 1 due to the kinematics of the non- wave-function vertex. We note that 
Eq. (||D essentially takes the same form as the LF bound-state equation (similar to the 
LF projection of BS equation) except the difference in kinematics(e.g. — x 2 > for the 



non- wave-function vertex). Incidentally, Einhorn |10| also discussed the extension of the 
LF BS amplitude in 1+1 QCD to a non-wave-function vertex similar to what we obtained 
in this work. 

In the above procedure, we also find that the four-body energy denominator D^ is ex- 
actly cancelled in the sum of LF time-ordered amplitudes as shown in Figs. 1(b) and 1(c), 
i.e., I/D4D2 + I/D4D2 = I/D2D2. We thus obtain the amplitude identical to the nonva- 
lence contribution in terms of ordinary LF wave functions of gauge boson(iy) and hadron 
(white blob) as drawn in Fig. 1(d). This method, however, requires to have some relevant 
operator depicted as the black square(/C) in Fig. 2(See also Fig. 1(d)), that is in general 
dependent on the involved momenta connecting one-body to three-body sector. We now 
present some details of kinematics in the semileptonic decay processes to discuss a rea- 
soning of how we handle the nonvalence contribution involving the momentum-dependent 
K, for relatively small momentum transfer processes such as 7r e 3, K^ and D — > Klv. 



In purely longitudinal momentum frame where q + > and Pu_ = P2_i_ = 0, the 
momentum transfer q 2 = q + q~~ can be written in terms of the momentum fraction a = 
P2 / P\ — 1 — Q + / P\ a s q 2 = (1 — a)(M 2 — Mf/a). With the iteration procedure in this 
frame, the results for the "+" -component of the current J^ are given by 
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/ ^r^y / ^^ k ^> y> u)*/(y, D, (5) 

where x = k + /P^, x' = x/a, D\ = q~ — V\ — (— pjj") with — p^ > and \I/jr/](x, kj_) = 
^i[2]/(^i 2 [2] - MW with M^jMo^] = m 2 ± /(l - x) + m 2 ± /x[m 2 ± /(l - x>) + m 2 ± /x'] 
and m 2 ± = m 2 + k\. Here, Eq. (0)(SD type equation) was folded in the derivation of 
Eq. flU) by the iteration procedure. While the relevant operator K, is in general dependent 
on all internal momenta (x, kj_,y, lj_), a sort of average on /C over y and lj_ in Eq.(5) 
which we define as Gp 1 p 2 = /[<iy|[<i 2 l_i_]/C(£, kj_; y, lj_)\l//(?/, lj_) is dependent only on x and 
k^. Now, the range of the momentum fraction x depends on the external momenta for 
the embedded states. As shown in Eq.(5), the lower bound of x for the kernel in the 
nonvalence contribution is given by a which has the value a = MijM\ at the maximum 
q 2 . As the mass difference between the primary and secondary mesons gets smaller, not 
only the range of q 2 is reduced but also a gets closer to 1. Perhaps, the best experimental 
process for such limit may be the pion beta decay ir^ — > ir°e ± u e , where our numerical 
prediction /_(0)//+(0) = —3.2 x 10~ 3 following the treatment presented in this work is 
in an excellent agreement with —3.5 x 10~ 3 obtained by the method proposed by Jaus 
including the zero-modes [11|. In Ademollo-Gatto's SU(3) limit f[2|| , the q 2 range of the 



nonvalence contribution shrinks to zero and a becomes precisely 1. However, even if a is 
not so close to 1, the initial wavefunction $?i(x, kjj plays the role of a weighting factor 
in the nonvalence contribution and enfeeble the contribution from the region of x near 
1. Thus, for the processes that we discussed in this letter, the effective x region for the 
nonvalence contribution is quite narrow. Similarly, the region of the transverse momentum 
kx is also limited only up to the scale of hadron size due to the same weighting factor 
\I/j(x, kj_). In this work, we thus approximate Gp 1 p 2 as a constant and examine the validity 
of this approximation by checking the frame independence of our numerical results. 



In Eqs.(4) and (5), the trace terms Sy(pg = k on ) = (4P 1 + /x / ){k^_ + [xm\ + (1 — 
x)m q ] [x'm 2 + (1 - x')m q ]} and S%y = S^(pi = p^ on ) + 4:pf on pt on (p^ - p^ on ) correspond to 
the product of initial and final LF spin-orbit wave functions that are uniquely determined 
by a generalized off-energy- shell Melosh transformation. Here, the subscript (on) means 
on-mass-shell and the instantaneous part of nonvalence diagram corresponds to S^ v — 
Sv(.P7 = Pion)- While the LF vertex function h^Z formally stems from H^X, practical 
informations on the radial wave function \I/j^](x, kj_) (consequently /if™) can be obtained 
by LF CQM. The details of our variational procedure to determine both mass spectra 
and wave functions of pseudoscalar mesons were recently documented in Refs. P,ffl along 
with an extensive test of the model in the spacelike exclusive processes. The same model 
is used in this work. 

For the check of frame-independence, we also compute the "+" component of the cur- 
rent Jjy in the Drell-Yan-West (q + = 0) frame where only valence contribution exists. 
Since the form factor f+(q 2 ) obtained from J^ in q + = frame is immune to the zero- 
mode contribution Pj7| JT3|JTll| , the comparison of f+(q 2 ) in the two completely different 
frames (i.e. q + = and q + ^ 0) would reveal the validity of existing model with respect 
to a covariance. The comparison of /_(g 2 ), however, cannot give a meaningful test of 



covariance because of the zero- mode complication as noted in Ref. [11]. Indeed, the dif- 
ference between the two (q + = and q + ^ 0) results of f-(q 2 ) amounts to the zero-mode 
contribution |13| . 



In our numerical calculation for the processes of K^ and D° -^ K l + v decays, we use 
the linear potential parameters presented in Ref. ||. In Fig. 3, we show the weak form 
factors f + (q 2 ) and f (q 2 ) = f+(q 2 ) +q 2 f-(q 2 )/(M 2 - M 2 ) for K$ 3 decays. The thick solid 
lines are our analytic solutions obtained from the q + = frame; note here again that 
the lower thick solid line (/o) in Fig. 3 is only the partial result without including the 
zero-mode contribution while the upper thick solid line (/ + immune to the zero-mode) is 
the full result. The thin solid lines are the full results of our effective calculations with a 
constant ((7^=3.95) fixed by the normalization of /+ at q 2 = limit. For comparison, 
we also show only the valence contributions(dotted lines) in q + ^ frame. As expected, 
a clearly distinguishable nonvalence contribution is found. Following the popular linear 
parametrization fH| , we plot the results of our effective solutions(thin solid lines) using 
fi(q 2 ) = fi(q 2 = w|)(l + ^i1 2 /M 2 + )(i = +,0). In comparison with the data, the same 
normalization as the data /+(0) = 1 |15j was used in Fig. 3. Our effective solution(upper 
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Figure 3. The weak form factors for K% compared with the experimental data [15 . 



thin solid line) is not only in a good agreement with the data [15] but also almost identical 
to that in q + = frame(upper thick solid line) indicating the frame-independence of our 
model. Note also that the difference in fo(q 2 ) between q + ^ (lower thin solid line) and 
q + = (lower thick solid line) frames amounts to the zero-mode contribution [|T3|1 . 

In comparison with experimental data, we summarized our results of several experi- 
mental observables in Table 1; i.e. the actual value of /+(0), the slopes A + [A ] of f+(q 2 ) 
[fo{q 2 )} at q 2 = 0, £i=/_(0)//+(0), and the decay rates T{K%) and r(K% 3 ). In the second 
column of Table 1, our full results including nonvalence contributions are presented along 
with the valence contributions in the square brackets. In the third column of Table 1, 
the results in q + = frame are presented with [without] the instantaneous part. As one 
can see in Table 1, adding the nonvalence contributions clearly improves the results of 
Ao, i.e. our full result of Ao= 0.025 is in an excellent agreement with the data, A xp ' = 
0.025±0.006. The results of £4 and F(K? l3 ) seem to be also improved by our full effective 
calculations including the nonvalence contributions. 

In Figs. 4(a,b), we show the weak form factors for D° — > K~£ + v decays and compare 
with the experimental data []l4|(full dot) with an error bar at Q 2 = as well as the lattice 



Table 1 

Model predictions for the parameters of Kg 3 decays. The decay width is in units of 10 6 
s _1 . The used CKM matrix is \V US \ = 0.2196 ± 0.0023 0. 





Effective 


q + = 


Experiment 


/+(0) 


0.962 [0.962] 


0.962 [0.962] 




A+ 


0.026 [0.083] 


0.026 [0.026] 


0.0288 ±0.0015[K e 3 ] 


Ao 


0.025 [-0.017] 


0.001 [-0.009] 


0.025 ± 0.006 [K°a] 


u 


-0.013 [-1.10] 


-0.29[-0.41] 


-0.11±0.09[/^ 3 ] 


r(4) 


7.3 ±0.15 


7.3 ±0.15 


7.5±0.08 


r(Ks) 


4.92 ±0.10 


4.66 ±0.10 


5.25±0.07 



QCD results []T6[ (circle and square) and |T7| (cross). All the line assignments are same as 
in Fig. 3. In Fig. 4(a), the thin solid line of our full result in q + ^ is not visible because 
it exactly coincides with the thick solid line of the result in q + = confirming the frame- 
independence of our calculations. Our value of /+(0)=0.736 is also within the error bar 
of the data |14[], / + xp '(0)=0.7±0.1. In Fig. 4(b), the difference between the thin and thick 
solid lines is the measure of the zero-mode contribution to fo(q 2 ) in q + = frame |]13|j . The 
form factors obtained from our effective calculations((j£>.ft;=3.5) are also plotted with the 
usual parametrization of pole dominance model, i.e. /+(o)(<7 2 ) = /+(o)(0)/(l — q 2 /Mf-^ Q+ ^). 
Our pole masses turn out to be Mi-=2.16 GeV and M +=2.79 GeV, respectively, and we 
note that Mi-=2.16 GeV is in a good agreement with the mass of -D*,i.e. 2.1 GeV. Using 
CKM matrix element \V CS \ = 1.04 ±0.16 0, our branching ratios Br(L>° 3 ) = 3.73 ± 1.24 



and Br(_D° 3 ) = 3.60 ± 1.19 are also comparable with the experimental data 3.64 ± 0.18 
and 3.22 ± 0.17 0, respectively. 

In summary, we presented an effective treatment of the LF nonvalence contributions 
crucial in the timelike exclusive processes. Using a SD-type approach and summing the LF 
time-ordered amplitudes, we obtained the nonvalence contributions in terms of ordinary 
LF wavefunctions of gauge boson and hadron that have been extensively tested in the 
spacelike exclusive processes. Including the nonvalence contribution, our results show a 
definite improvement in comparison with experimental data on K^ and D° — > K~£ + i>e 
decays. Our result on 7r e3 is also consistent with the result obtained by other method. 
Furthermore, the frame-independence of our results indicate that a constant Gp 1 p 2 is an 
approximation appropriate to the small momentum transfer processes. Applications to 
the heavy-to-light decay processes involving large momentum transfers would require an 
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Figure 4. The weak form factors for D — > K transition. 
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improvement on this approximation perhaps guided by the perturbative QCD approach. 
Consideration along this line is underway. 
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